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We present large-scale Monte Carlo results for the dynamical critical exponent z and the spatio- 
temporal two-point correlation function of a (2+l)-dimensional quantum XY model with bond 
dissipation, proposed to describe a quantum critical point in high-T c cuprates near optimal doping. 
The phase variables of the model, originating with a parametrization of circulating currents within 
the CUO2 unit cells in cuprates, are compact, {0 r>r } G [— tt,tt). The dynamical critical exponent is 
found to be z ~ 1, and the spatio-temporal correlation functions are explicitly demonstrated to be 
isotropic in space-imaginary time. The model thus has a fluctuation spectrum where momentum 
and frequency enter on equal footing, rather than having the essentially momentum-independent 
marginal Fermi liquid-like fluctuation spectrum previously reported for the same model. 
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Quantum critical points describe systems with diverg- 
ing length scales at zero temperature, and have come 
into much focus in recent years as possible descriptions 
of anomalous phenomena in strongly correlated fermion 
systems and systems with competing orders. 1 One prime 
example of this is represented by the high-T c supercon- 
ducting cuprates, where various types of quantum critical 
phenomena have been proposed as essential for under- 
standing the many unusual normal-state transport prop- 
erties these systems exhibit. This has, over the past quar- 
ter of a century, represented one of the major challenges 
in condensed-matter physicsP 

One successful phenomenological framework is to de- 
scribe the the normal phase around optimal doping as 
a marginal Fermi liquid (MFL), 3 the weakest possible 
violation of having a non-zero quasiparticle residue at 
the Fermi surface. Among the merits of the MFL phe- 
nomenology is that it describes transport properties in 
this strange metallic phase in good accordance with ex- 
periments. This follows naturally from the essentially 
momentum-independent, linear-in- frequency, fluctuation 
spectrum of the MFL hypothesisP 

More recent works have pursued a more microscopic 
foundation of MFL. The underlying picture is that there 
exists a quantum critical point (QCP) residing at T = 
beneath the superconducting dome 3 The degrees of free- 
dom associated with this QCP are circulating currents 
within the unit cells of the Cu02 layers. The main 
idea is that the MFL phenomenology arises from the 
quantum critical fluctuations of these currents above 
the QCP at T > T c . It has also been demonstrated 
how the same fluctuations may give rise to e?-wave high- 
Tc superconductivity.^ The ordering of such circulat- 
ing currents upon lowering the temperature from the 
strange metal region into the pseudogap region is a can- 
didate for a possible competing order in this part of the 
phase diagram. - Magnetic order conforming with such 
circulating currents has in fact been observed in sev- 
eral experiments.^^ It must be mentioned that others 
argue that such signatures may have a quite different 



originJ^EU and also numerical results disagree on the 
presence of such circulating currents] EMU] b u t the model 
remains one of t he central theories of the physics of high- 
Tc cuprates.^ 

A remarkable implication of a (/-independent fluctua- 
tion spectrum, such as that posited in MFL theory, is 
that the associated QCP exhibits local quantum criti- 
cality (LQC). Defining the dynamical critical exponent 
z from the scaling of momentum and frequency at the 
quantum critical point, uj ~ q z , this means that, for- 
mally, z = oo. It is a highly nontrivial question as to 
how such a remarkable property of a quantum critical 
point can a rise in an extended system. Recently, it was 
argue d 20 * 21 ^ that precisely such local criticality is found 
in a (2+l)-dimensional quantum XY model with bond 
dissipation of the Caldeira-LeggettP^ form. The angle 
variables of this model were associated with circulating 
current degrees of freedom, as will be explained below. 

The results of Ref . 1201 would imply that the previ- 
ously hypothesized MFL fluctuation spectrum has been 
derived from a microscopic theory applicable to cuprates. 
In a broader perspective, it is of considerable interest to 
investigate in detail if such unusual behavior can occur in 
model systems of condensed matter, as related variants 
of locality have also been considered in the context of 
gauge/gravity duality 23 and QCPs in disordered systems 
and heavy fermion compounds.^ 

From naive scaling argument d 24 * 25 ! applied to the dis- 
sipative model proposed in Ref. [2111 one might expect 
that dissipation is irrelevant in the renormalization group 
sense. The result would then not be LQC, but instead 
conventional quantum criticality with isotropic scaling, 
z = 1. Here, we report results from Monte Carlo simula- 
tions performed directly on the (2+l)-dimcnsional quan- 
tum XY model with bond dissipation and compact angle 
variables, considered in Ref. [20] Our results strongly in- 
dicate that in this model, z = 1. 

The dissipative (2+l)-dimensional [(2+l)D] XY ac- 
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tion considered in Ref. 20J takes the form 

S = -K J2 C° S (A6» r ,r',r) (1) 

<r,r') t=1 

- Kt ^2 ^2 COs(6> r , r+ i - r , T ) 
r t=1 

a sr-^ ( ir \ 2 (A6* r r ' !r — A6*r jr ' iT ') 2 

+ 2~At", VW sin 2 (^|r~r'|) ' 

<r,r') t^t' x lt i 

when put on a cubic L x L x L T lattice. The bond vari- 
ables are given by A# r . r ' jT = 6 r T — 9 r ' tT , where the sum 
over r and r' goes over nearest neighbors in the x-y plane. 
Periodic boundary conditions arc implicit in the imagi- 
nary time direction, and are also applied in the spatial 
directions. 

Such a model has previously been employed as an effec- 
tive description of a resistively shunted Josephson junc- 
tion arrayJ^Sl and it may also be viewed as a generic quan- 
tum rotor model with dissipative currents. A third pos- 
sible interpretation in the context of high-T c cuprates 
is as follows. Suppose the angles a priori can take 
only four possible values. These four values then repre- 
sent the directions of a pseudospin associated with the 
four possible ordered circulating current patterns within 
each Cu0 2 unit cell (see, e.g., Fig. 1 of Ref. [27J. The 
first two terms represent the standard interaction ener- 
gies in space-imaginary time of these circulating currents 
in neighboring unit cells, and have been derived from 
microscopies The last term is the term responsible for 
dissipating the ordered circulating currentsPEl 

In Eq. ([I]), the angles are continuous variables. We 
will discuss a possible a posteriori justification for this 
later in this Rapid Communication, by showing that an 
added four-fold anisotropy term is perturbatively irrele- 
vant. Reference HOJ moreover, appears to treat r ^ T as 
compact variables, also in the presence of a dissipation 
term that apparently renders the action nonperiodic in 
the angle variables! 2 ^ In order to investigate numerically 
the same model considered in Ref. [501 we therefore 
compactify the expression A6 r y_ T — A6 r y_ T i so that it 
is defined modulo 2ir. We will discuss alternative choices 
later. 

The calculations of Ref. [20] were not restricted to 
any specific parameter regime, but predicted that every 
point on the T = quantum critical surface in a—K—K T 
(parameter) space (for a > 0) should be a local QCP. Ac- 
cordingly, we choose convenient coupling constants when 
searching for LQC in our simulations, and for the re- 
sults presented here, the dissipation strength is fixed at 
a = 0.05. 

The phase diagram (not shown) is qualitatively very 
similar to those found for relat ed compact (1+1)D mod- 
els with bond dissipation! 25 * 29 ! It features a single criti- 
cal surface that separates a disordered from a fully or- 
dered phase, and which is continuously connected to the 
3D XY critical line at a = 0. For similar models in 



(1+1) dimensions, only the region of relatively moder- 
ate dissipation was accessible to simulations, as increas- 
ing a increases finite-size effects, resulting in apparent 
values z < 1 for the dynamical critical exponent. As 
expected, this problem is no less severe in (2+1) dimen- 
sions. Available system sizes are restricted by the absence 
of cluster algorithms to treat models with bond dissi- 
pation appropriately} 2 ^ and we are therefore confined to 
local Metropolis updates. 

To locate the phase transition, we vary the spatial cou- 
pling K and use the crossing point for different system 
sizes L of the Binder cumulant g = 1— (|m| 4 )/(2(|m| 2 ) 2 ). 
Here, m = J2r r ex P t^r.r] is the order parameter of 
the t/(l)-symmetric degrees of freedom. Due to the 
anisotropy of the interactions, we have to calculate g for 
multiple values of L T for each spatial system size L, as 
described in more detail, e.g., in Ref. 25] The value 
L T = L* T where the function g(L T ) reaches its maximum 
corresponds to the optimal temporal extent for which the 
system appears as isotropic as it can be, the anisotropic 
interactions taken into account. 

For a conventional QCP, at which the correlation 
length £ r in imaginary time scales with the correlation 
length £ in space as £ T ~ £ z with a finite z, we expect 
to observe the scaling relation L* ~ L z . This scaling 
procedure then allows one to extract the dynamical crit- 
ical exponent z from Binder cumulant data. For a local 
QCP formally having z — oo, we expect this scaling to 
break down. Our strategy to search for possible LQC in 
the model (JT|) is therefore to perform the above proce- 
dure assuming conventional criticality, and then look for 
indications that this hypothesis should be rejected. 
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FIG. 1: (Color online) Finite-size analysis of the maximum 
L* of the Binder cumulant curves g(L T ) as a function of spa- 
tial system size L. For the black data points, the dynamical 
critical exponents z as given in Table [I] are obtained from the 
slope of the fitting lines (dashed). The red (gray) points show 
similar results for site dissipation for comparison, where a fit 
of the three largest systems yields z = 1.84(3). 

The results of this finite-size analysis is shown in Fig. 
[T] with the values of the dynamical critical exponent z 
given in Table [T] Here, we have chosen three different 
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TABLE I: Critical coupling K c and dynamical critical expo- 
nent z for different values of the quantum coupling K T , but 
for the same dissipation strength a = 0.05. Uncertainty es- 
timates for z have been calculated by a bootstrap procedure, 
including the uncertainty in K c . 





K c 


z 


0.2 


0.48068(5) 


0.968(8) 


0.6 


0.28244(4) 


0.985(8) 


1.0 


0.18008(5) 


0.970(11) 



values of the quantum coupling K T in order to investi- 
gate both the limit of relatively weak quantum coupling 
and the opposite limit leading to relatively strong system 
anisotropy. 

The results show that the effective dynamical critical 
exponent is z < 1 for all the parameter sets considered, 
and we expect that we could obtain z w 1 if we were 
able to reach higher values of L. (For a smaller value 
a = 0.02, we obtained z = 1 within statistical uncer- 
tainty.) It is conceivable that signatures of LQC would 
be visible only for systems larger than the admittedly 
moderate system sizes accessible to present algorithms. 
However, were that the case, the true z — > oo nature of 
the model would likely reveal itself as strongly increasing 
effective values of z as a finite-size effect for increasing 
L. For comparison, we have also carried out simulations 
with equivalent parameters of a (2+l)D XY model with 
site dissipation, for which z = 2 is expected! 25 * 30 ! The 
results are included in Fig. [T] and already for system 
sizes comparable to those for bond dissipation, we ob- 
serve (finite-size) crossover behavior with z — > 2. For 
bond dissipation, we observe no tendency toward z > 1 
for either of the parameter sets, and it is hard to imagine 
how crossover to z — > oo scaling should be much slower 
than crossover to z — 2 scaling. 

For all results reported here, we have used parallel 
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tempering*^ to reduce autocorrelation times, and to en- 
sure that the simulations are well equilibrated. To emu- 
late the continuous U(l) symmetry, the simulations are 
made for Z q clock models, with q = 128 for K T = 0.2, 1.0, 
and q = 32 for K T = 0.6. The nature of the criticality 
remains unchanged also when increasing to q = 1024. 
The results are obtained using an implementation of the 
Mersenne Twister^ random number generator, but other 
random number generators produced consistent results. 

Although we found no indication of LQC from the scal- 
ing of the Binder cumulant, we also considered the cor- 
relations of the order parameter field directly, 



C(r 



r') 



< T e~ 



(2) 



The correlation functions presented here are obtained for 
the parameter set K T = 0.6, with L T = L* and K = K c 
as obtained from the previous simulations, and therefore 
serve as a self-consistency check of the Binder scaling 



procedure. From Fig. [2] it is evident that the correlation 
function at the critical point decays isotropically in space- 
imaginary time. In other words, there are no signs of 
locality. 

Furthermore, we have verified that the same conclusion 
may be drawn for the other values of K T considered, and 
also for larger system sizes with aspect ratios found from 
extrapolation based on the power law shown in Fig. [T] 
As an additional test, we compared the correlation func- 
tions shown here with those obtained by setting a = 
in Eq. ([lj. Letting K T > K, values of L T and K c were 
determined by the same procedure as for the dissipative 
model. There is no indication that adding dissipation 
changes the scaling of the temporal correlation length £ r 
with respect to the spatial correlation length £. 
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FIG. 2: (Color online) Correlation functions at the critical 
point K c = 0.28244 for dissipation strength a = 0.05 and 
quantum coupling K T = 0.6. The system size L = 32, 
L T = 49 as L* corresponds to the rightmost data point of the 
midmost data series in Fig. [T] The correlation function is de- 
fined in the spatial direction as g(x) = g(|r— r'|) = C(r— r',0) 
and in the temporal direction as g(j) = C(0,r), with C de- 
fined in Eq. |5J. Also, L x = L. Error bars are smaller than 
the line width, and the dotted lines are guides to the eye. 

Depending on how A0 r r /. T is interpreted in the dis- 
sipation term, it may be argued either that the correct 
treatment is to compactify only the gradients A9 r r r tT , re- 
stricting them to the interval [— n, 7r), or to do so to the 
difference A0 r-r ' iT — A0 T ^iy as well. Although we have 
chosen the latter, as in Ref. HHJ we also performed simu- 
lations with the former compactification scheme. The re- 
sults are qualitatively similar, with the difference merely 
amounting to a renormalization of the dissipative cou- 
pling a. In other words, the absence of LQC in this 
model is not contingent on the choice of compactification 
scheme. 

As explained in connection with Eq. (JT|), the underly- 
ing circulating current degrees of freedom are most nat- 
urally described by discrete, .^-symmetric variables. In 
Ref. 21, it was argued that a model with continuous 
C/(l) symmetry nonetheless would be a correct descrip- 
tion. The result of LQC would then also apply to the 
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four-state model of the original degrees of freedom since 
a fourfold anisotropy field, given by 



S4 = hi cos (40 r)T ), 



(3) 



would be irrelevant at the critical point of the action (JT|) . 
We have investigated the effect of a fourfold anisotropy in 
our simulations by including the term (3| in the action. 
Using the approach of Ref. ESI we find the same result 
for the dissipative (2+l)D XY model as reported there 
for the classical 3D XY model, namely, that the /14 term 
is perturbatively irrelevant. 

The soft constraint represented by a (finite) anisotropy 
term is not obviously the same as the hard constraint con- 
stituted by the discrete Z4 variables of the original model 
(the limit h± = 00). We may only speculate whether a 
putative LQC fixed point for a U{1) theory might sur- 
vive in the limit /14 —> 00, but note that our simulations 
showed no signs of locality neither when enforcing a soft 
nor a strong Z4 constraint on the variables.^ 

Finally, we briefly consider variants of LQC other than 
that of Ref. HTJ which predicts a strictly infinite z for 
£ r ~ £ z so that £ is strictly vanishing at criticality. An- 
other conceivable sense in which z — > 00 is by activated 
dynamical scaling,^ i.e., scaling on the form ln£ r ~ 



In this case, as we expect also in the first case, local- 
ity would manifest itself in our simulations as a strongly 
increasing value of z > 1 as the thermodynamical limit 
was approached. This is not observed in our results. We 
have also verified explicitly, by an appropriate modifica- 
tion of the scalingpSl that our results are not consistent 
with activated dynamical scaling. 

In conclusion, we find no signs of local quantum crit- 
icality in the compact (2+l)D XY model with bond 
dissipation, but instead conventional quantum criticality 
with indications of isotropic scaling of imaginary time 
and space. This implies that the fluctuation spectrum 
of the model is a function of the combination q 2 + u 2 , 
rather than being dependent only on the frequency uj, 
but not on the momentum q (which would be a hallmark 
of local quantum criticality) . Our results therefore differ 
in a fundamental way from those obtained from the same 
model in Ref. [201 
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